Transcriptomic and functional genetic evidence for distinct ecophysiological responses across complex life cycle stages.

Organisms with complex life cycles demonstrate a remarkable ability to change their phenotypes across development, presumably as an evolutionary adaptation to developmentally variable environments. Developmental variation in environmentally sensitive performance, and thermal sensitivity in particular, has been well documented in holometabolous insects. For example, thermal performance in adults and juvenile stages exhibit little genetic correlation (genetic decoupling) and can evolve independently, resulting in divergent thermal responses. Yet, we understand very little about how this genetic decoupling occurs. We tested the hypothesis that genetic decoupling of thermal physiology is driven by fundamental differences in physiology between life stages, despite a potentially conserved cellular stress response. We used RNAseq to compare transcript expression in response to a cold stressor in Drosophila melanogaster larvae and adults and used RNA interference (RNAi) to test whether knocking down nine target genes differentially affected larval and adult cold tolerance. Transcriptomic responses of whole larvae and adults during and following exposure to -5°C were largely unique both in identity of responding transcripts and in temporal dynamics. Further, we analyzed the tissue-specificity of differentially expressed transcripts from FlyAtlas 2 data, and concluded that stage-specific differences in transcription were not simply driven by differences in tissue composition. In addition, RNAi of target genes resulted in largely stage-specific and sometimes sex-specific effects on cold tolerance. The combined evidence suggests that thermal physiology is largely stage-specific at the level of gene expression, and thus natural selection may be acting on different loci during the independent thermal adaptation of different life stages.

Stage-specific genotype-by-environment interactions for cold and heat hardiness in Drosophila melanogaster.

Environments often vary across a life cycle, imposing fluctuating natural selection across development. Such fluctuating selection can favor different phenotypes in different life stages, but stage-specific evolutionary responses will depend on genetic variance, covariance, and their interaction across development and across environments. Thus, quantifying how genetic architecture varies with plastic responses to the environment and across development is vital to predict whether stage-specific adaptation will occur in nature. Additionally, the interaction of genetic variation and environmental plasticity (GxE) may be stage-specific, leading to a three-way interaction between genotype, environment, and development or GxDxE. To test for these patterns, we exposed larvae and adults of Drosophila melanogaster isogenic lines derived from a natural population to extreme heat and cold stress after developmental acclimation to cool (18 °C) and warm (25 °C) conditions and measured genetic variance for thermal hardiness. We detected significant GxE that was specific to larvae and adults for cold and heat hardiness (GxDxE), but no significant genetic correlation across development for either trait at either acclimation temperature. However, cross-development phenotypic correlations for acclimation responses suggest that plasticity itself may be developmentally constrained, though rigorously testing this hypothesis requires more experimentation. These results illustrate the potential for stage-specific adaptation within a complex life cycle and demonstrate the importance of measuring traits at appropriate developmental stages and environmental conditions when predicting evolutionary responses to changing climates.

Phenology of Drosophila species across a temperate growing season and implications for behavior.

Drosophila community composition is complex in temperate regions with different abundance of flies and species across the growing season. Monitoring Drosophila populations provides insights into the phenology of both native and invasive species. Over a single growing season, we collected Drosophila at regular intervals and determined the number of individuals of the nine species we found in Kansas, USA. Species varied in their presence and abundance through the growing season with peak diversity occurring after the highest seasonal temperatures. We developed models for the abundance of the most common species, Drosophila melanogaster, D. simulans, D. algonquin, and the recent invasive species, D. suzukii. These models revealed that temperature played the largest role in abundance of each species across the season. For the two most commonly studied species, D. melanogaster and D. simulans, the best models indicate shifted thermal optima compared to laboratory studies, implying that fluctuating temperature may play a greater role in the physiology and ecology of these insects than indicated by laboratory studies, and should be considered in global climate change studies.

Speed of exposure to rapid cold hardening and genotype drive the level of acclimation response in Drosophila melanogaster.

Comparing methodologies that attempt to mimic natural conditions is important when evaluating thermal tolerances of ectotherms, as exposing animals to different artificial thermal regimes may provide conflicting information of an insect's thermal profile. Rapid cold hardening (RCH) occurs in ectotherms and typically increases survivorship to extreme cold exposure through a short, pre-treatment to a non-lethal cold temperature. Here we assess survivorship in a set of genotypes from the Drosophila melanogaster Reference Panel for direct and ramping RCH pre-treatments at cooling rates occurring under more natural conditions (0.1 °C/min and 0.5 °C/min) in combination with a direct and ramped rewarming treatment post cold exposure. We find that all three pre-treatment exposures alone significantly increase survivorship. We find significant correlations in survivorship among treatments across genotypes, suggesting that regardless of the pre-treatment, individuals of a given genotype have an innate level of acclimation. When rewarming is introduced, survivorship significantly decreased relative to pre-treatment alone and correlations of survival between phenotypes were not significant. Our results suggest that rewarming and slow RCH are costly to survival while a quicker RCH may impart physiological benefits more consistently across genotypes.

Costs of cold acclimation on survival and reproductive behavior in Drosophila melanogaster.

Fitness is determined by the ability of an organism to both survive and reproduce; however, the mechanisms that lead to increased survival may not have the same effect on reproductive success. We used nineteen natural Drosophila melanogaster genotypes from the Drosophila Genetic Reference Panel to determine if adaptive plasticity following short-term acclimation through rapid cold-hardening (RCH) affects mating behavior and mating success. We confirmed that exposure to the acclimation temperature is beneficial to survival following cold stress; however, we found that this same acclimation temperature exposure led to less efficient male courtship and a significant decrease in the likelihood of mating. Cold tolerance and the capacity to respond plastically to cold stress were not correlated with mating behavior following acclimation, suggesting that the genetic control of the physiological effects of the cold temperature exposure likely differ between survival and behavioral responses. We also tested whether the exposure of males to the acclimation temperature influenced courtship song. This exposure again significantly increased courtship duration; however, courtship song was unchanged. These results illustrate costs of short-term acclimation on survival and reproductive components of fitness and demonstrate the pronounced effect that short-term thermal environment shifts can have on reproductive success.

Cold adaptation does not alter ATP homeostasis during cold exposure in Drosophila melanogaster.

In insects and other ectotherms, cold temperatures cause a coma resulting from loss of neuromuscular function, during which ionic and metabolic homeostasis are progressively lost. Cold adaptation improves homeostasis during cold exposure, but the ultimate targets of selection are still an open question. Cold acclimation and adaptation remodels mitochondrial metabolism in insects, suggesting that aerobic energy production during cold exposure could be a target of selection. Here, we test the hypothesis that cold adaptation improves the ability to maintain rates of aerobic energy production during cold exposure by using 31 P NMR on live flies. Using lines of Drosophila melanogaster artificially selected for fast and slow recovery from a cold coma, we show that cold exposure does not lower ATP levels and that cold adaptation does not alter aerobic ATP production during cold exposure. Cold-hardy and cold-susceptible lines both experienced a brief transition to anaerobic metabolism during cooling, but this was rapidly reversed during cold exposure, suggesting that oxidative phosphorylation was sufficient to meet energy demands below the critical thermal minimum, even in cold-susceptible flies. We thus reject the hypothesis that performance under mild low temperatures is set by aerobic ATP supply limitations in D. melanogaster, excluding oxygen and capacity limitation as a weak link in energy supply. This work suggests that the modulations to mitochondrial metabolism resulting from cold acclimation or adaptation may arise from selection on a biosynthetic product(s) of those pathways rather than selection on ATP supply during cold exposure.

Antagonistic pleiotropy and mutation accumulation contribute to age-related decline in stress response.

As organisms age, the effectiveness of natural selection weakens, leading to age-related decline in fitness-related traits. The evolution of age-related changes associated with senescence is likely influenced by mutation accumulation (MA) and antagonistic pleiotropy (AP). MA predicts that age-related decline in fitness components is driven by age-specific sets of alleles, nonnegative genetic correlations within trait across age, and an increase in the coefficient of genetic variance. AP predicts that age-related decline in a trait is driven by alleles with positive effects on fitness in young individuals and negative effects in old individuals, and is expected to lead to negative genetic correlations within traits across age. We build on these predictions using an association mapping approach to investigate the change in additive effects of SNPs across age and among traits for multiple stress-response fitness-related traits, including cold stress with and without acclimation and starvation resistance. We found support for both MA and AP theories of aging in the age-related decline in stress tolerance. Our study demonstrates that the evolution of age-related decline in stress tolerance is driven by a combination of alleles that have age-specific additive effects, consistent with MA, as well as nonindependent and antagonistic genetic architectures characteristic of AP.

Genetic Decoupling of Thermal Hardiness across Metamorphosis in Drosophila melanogaster.

As organisms age the environment fluctuates, exerting differential selection across ontogeny. In particular, highly seasonal environments expose life stages to often drastically different thermal environments. This developmental variation is particularly striking in organisms with complex life cycles, wherein life history stages also exhibit distinct morphologies, physiologies, and behaviors. Genes acting pleiotropically on thermal responses may produce genetic correlations across ontogeny, constraining the independent evolution of each life stage to their respective thermal environments. To investigate whether developmental genetic correlations constrain the evolution thermal hardiness of the fly Drosophila melanogaster, we applied quantitative genetic analyses to cold hardiness measured in both larvae and adults from isogenic lines of the Drosophila Genetic Reference Panel (DGRP), using survival at stressful low temperatures as the phenotypic metric. Using full genome resequencing data for the DGRP, we also implemented genome-wide association (GWA) analysis using Bayesian Sparse Linear Mixed Models (BSLMMs) to estimate associations between naturally segregating variation and cold hardiness for both larvae and adults. Quantitative genetic analyses revealed no significant genetic correlation for cold hardiness between life stages, suggesting complete genetic decoupling of thermal hardiness across the metamorphic boundary. Both quantitative genetic and GWA analyses suggested that polygenic variation underlies cold hardiness in both stages, and that associated loci largely affected one stage or the other, but not both. However, reciprocal enrichment tests and correlations between BSLMM parameters for each life stage support some shared physiological mechanisms that may reflect common cellular thermal response pathways. Overall, these results suggest no developmental genetic constraints on cold hardiness across metamorphosis in D. melanogaster, an important consideration in evolutionary models of responses to changing climates. Genetic correlations for environmental sensitivity across ontogeny remains largely unexplored in other organisms, thus assessing the generality of genetic decoupling will require further quantitative or population genetic analysis in additional species.

Seasonal variation in basal and plastic cold tolerance: Adaptation is influenced by both long- and short-term phenotypic plasticity.

Understanding how thermal selection affects phenotypic distributions across different time scales will allow us to predict the effect of climate change on the fitness of ectotherms. We tested how seasonal temperature variation affects basal levels of cold tolerance and two types of phenotypic plasticity in Drosophila melanogaster. Developmental acclimation occurs as developmental stages of an organism are exposed to seasonal changes in temperature and its effect is irreversible, while reversible short-term acclimation occurs daily in response to diurnal changes in temperature. We collected wild flies from a temperate population across seasons and measured two cold tolerance metrics (chill-coma recovery and cold stress survival) and their responses to developmental and short-term acclimation. Chill-coma recovery responded to seasonal shifts in temperature, and phenotypic plasticity following both short-term and developmental acclimation improved cold tolerance. This improvement indicated that both types of plasticity are adaptive, and that plasticity can compensate for genetic variation in basal cold tolerance during warmer parts of the season when flies tend to be less cold tolerant. We also observed a significantly stronger trade-off between basal cold tolerance and short-term acclimation during warmer months. For the longer-term developmental acclimation, a trade-off persisted regardless of season. A relationship between the two types of plasticity may provide additional insight into why some measures of thermal tolerance are more sensitive to seasonal variation than others.

Age-Related Changes in Locomotor Performance Reveal a Similar Pattern for Caenorhabditis elegans, Mus domesticus, Canis familiaris, Equus caballus, and Homo sapiens.

Locomotion is one of the major physiological functions for most animals. Previous studies have described aging mechanisms linked to locomotor performance among different species. However, the precise dynamics of these age-related changes, and their interactions with development and senescence, are largely unknown. Here, we use the same conceptual framework to describe locomotor performances in Caenorhabditis elegans, Mus domesticus, Canis familiaris, Equus caballus, and Homo sapiens. We show that locomotion is a consistent biomarker of age-related changes, with an asymmetrical pattern throughout life, regardless of the type of effort or its duration. However, there is variation (i) among species for the same mode of locomotion, (ii) within species for different modes of locomotion, and (iii) among individuals of the same species for the same mode of locomotion. Age-related patterns are modulated by genetic (such as selective breeding) as well as environmental conditions (such as temperature). However, in all cases, the intersection of the rising developmental phase and the declining senescent phase reveals neither a sharp transition nor a plateau, but a smooth transition, emphasizing a crucial moment: the age at peak performance. This transition may define a specific target for future investigations on the dynamics of such biological interactions.

Cold adaptation increases rates of nutrient flow and metabolic plasticity during cold exposure in Drosophila melanogaster.

Metabolic flexibility is an important component of adaptation to stressful environments, including thermal stress and latitudinal adaptation. A long history of population genetic studies suggest that selection on core metabolic enzymes may shape life histories by altering metabolic flux. However, the direct relationship between selection on thermal stress hardiness and metabolic flux has not previously been tested. We investigated flexibility of nutrient catabolism during cold stress in Drosophila melanogaster artificially selected for fast or slow recovery from chill coma (i.e. cold-hardy or -susceptible), specifically testing the hypothesis that stress adaptation increases metabolic turnover. Using (13)C-labelled glucose, we first showed that cold-hardy flies more rapidly incorporate ingested carbon into amino acids and newly synthesized glucose, permitting rapid synthesis of proline, a compound shown elsewhere to improve survival of cold stress. Second, using glucose and leucine tracers we showed that cold-hardy flies had higher oxidation rates than cold-susceptible flies before cold exposure, similar oxidation rates during cold exposure, and returned to higher oxidation rates during recovery. Additionally, cold-hardy flies transferred compounds among body pools more rapidly during cold exposure and recovery. Increased metabolic turnover may allow cold-adapted flies to better prepare for, resist and repair/tolerate cold damage. This work illustrates for the first time differences in nutrient fluxes associated with cold adaptation, suggesting that metabolic costs associated with cold hardiness could invoke resource-based trade-offs that shape life histories.

Artificial selection on chill-coma recovery time in Drosophila melanogaster: Direct and correlated responses to selection.

Artificial selection can be used to create populations with extreme phenotypic responses to environmental stressors. When artificial selection is applied to a single component of a stress response, this selection may result in correlated responses in other stress responses, a phenomenon called cross-tolerance, which is ultimately controlled by the genetic correlations among traits. We selected for extreme responses to cold tolerance by selecting for chill-coma recovery time from a single temperate population of Drosophila melanogaster. Chill-coma recovery time is a common metric of low, but non-lethal, cold temperature tolerance. Replicated divergent artificial selection was applied to a genetically variable base population for 31 generations, resulting in two cold resistant, two cold susceptible, and two unselected control lines. To quantify the relationship between selection on chill-coma recovery and other metrics of thermal performance, we also measured survivorship after acute cold exposure, survivorship after chronic cold exposure, survivorship after cold exposure following a pre-treatment period (rapid cold hardening), starvation tolerance, and heat tolerance. We find that chill-coma recovery time is heritable within this population and that there is an asymmetric response to increased and decreased chill-coma recovery time. Surprisingly, we found no cross-tolerances between selection on chill-coma recovery time and the other environmental stress response traits. These results suggest that although artificial selection has dramatically altered chill-coma recovery time, the correlated response to selection on other stress response phenotypes has been negligible. The lack of a correlated response suggests that chill-coma recovery time in these selection lines is likely genetically independent from measures of cold survivorship tested here.

Adaptation to Low Temperature Exposure Increases Metabolic Rates Independently of Growth Rates.

Metabolic cold adaptation is a pattern where ectotherms from cold, high-latitude, or -altitude habitats have higher metabolic rates than ectotherms from warmer habitats. When found, metabolic cold adaptation is often attributed to countergradient selection, wherein short, cool growing seasons select for a compensatory increase in growth rates and development times of ectotherms. Yet, ectotherms in high-latitude and -altitude environments face many challenges in addition to thermal and time constraints on lifecycles. In addition to short, cool growing seasons, high-latitude and - altitude environments are characterized by regular exposure to extreme low temperatures, which cause ectotherms to enter a transient state of immobility termed chill coma. The ability to resume activity quickly after chill coma increases with latitude and altitude in patterns consistent with local adaptation to cold conditions. We show that artificial selection for fast and slow chill coma recovery among lines of the fly Drosophila melanogaster also affects rates of respiratory metabolism. Cold-hardy fly lines, with fast recovery from chill coma, had higher respiratory metabolic rates than control lines, with cold-susceptible slow-recovering lines having the lowest metabolic rates. Fast chill coma recovery was also associated with higher respiratory metabolism in a set of lines derived from a natural population. Although their metabolic rates were higher than control lines, fast-recovering cold-hardy lines did not have faster growth rates or development times than control lines. This suggests that raised metabolic rates in high-latitude and -altitude species may be driven by adaptation to extreme low temperatures, illustrating the importance of moving "Beyond the Mean".

Intraspecific variation of a dominant grass and local adaptation in reciprocal garden communities along a US Great Plains' precipitation gradient: implications for grassland restoration with climate change.

Identifying suitable genetic stock for restoration often employs a 'best guess' approach. Without adaptive variation studies, restoration may be misguided. We test the extent to which climate in central US grasslands exerts selection pressure on a foundation grass big bluestem (Andropogon gerardii), widely used in restorations, and resulting in local adaptation. We seeded three regional ecotypes of A. gerardii in reciprocal transplant garden communities across 1150 km precipitation gradient. We measured ecological responses over several timescales (instantaneous gas exchange, medium-term chlorophyll absorbance, and long-term responses of establishment and cover) in response to climate and biotic factors and tested if ecotypes could expand range. The ecotype from the driest region exhibited greatest cover under low rainfall, suggesting local adaptation under abiotic stress. Unexpectedly, no evidence for cover differences between ecotypes exists at mesic sites where establishment and cover of all ecotypes were low, perhaps due to strong biotic pressures. Expression of adaptive differences is strongly environment specific. Given observed adaptive variation, the most conservative restoration strategy would be to plant the local ecotype, especially in drier locations. With superior performance of the most xeric ecotype under dry conditions and predicted drought, this ecotype may migrate eastward, naturally or with assistance in restorations.

Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long- and short-term thermal acclimation.

Seasonal and daily thermal variation can limit species distributions because of physiological tolerances. Low temperatures are particularly challenging for ectotherms, which use both basal thermotolerance and acclimation, an adaptive plastic response, to mitigate thermal stress. Both basal thermotolerance and acclimation are thought to be important for local adaptation and persistence in the face of climate change. However, the evolutionary independence of basal and plastic tolerances remains unclear. Acclimation can occur over longer (seasonal) or shorter (hours to days) time scales, and the degree of mechanistic overlap is unresolved. Using a midlatitude population of Drosophila melanogaster, we show substantial heritable variation in both short- and long-term acclimation. Rapid cold hardening (short-term plasticity) and developmental acclimation (long-term plasticity) are positively correlated, suggesting shared mechanisms. However, there are independent components of these traits, because developmentally acclimated flies respond positively to short-term acclimation. A strong negative correlation between basal cold tolerance and developmental acclimation suggests that basal cold tolerance may constrain developmental acclimation, whereas a weaker negative correlation between basal cold tolerance and short-term acclimation suggests less constraint. Using genome-wide association mapping, we show the genetic architecture of rapid cold hardening and developmental acclimation responses are nonoverlapping at the SNP and corresponding gene level. However, genes associated with each trait share functional similarities, including genes involved in apoptosis and autophagy, cytoskeletal and membrane structural components, and ion binding and transport. These results indicate substantial opportunity for short-term and long-term acclimation responses to evolve separately from each other and for short-term acclimation to evolve separately from basal thermotolerance.

Ecotypes of an ecologically dominant prairie grass (Andropogon gerardii) exhibit genetic divergence across the U.S. Midwest grasslands' environmental gradient.

Big bluestem (Andropogon gerardii) is an ecologically dominant grass with wide distribution across the environmental gradient of U.S. Midwest grasslands. This system offers an ideal natural laboratory to study population divergence and adaptation in spatially varying climates. Objectives were to: (i) characterize neutral genetic diversity and structure within and among three regional ecotypes derived from 11 prairies across the U.S. Midwest environmental gradient, (ii) distinguish between the relative roles of isolation by distance (IBD) vs. isolation by environment (IBE) on ecotype divergence, (iii) identify outlier loci under selection and (iv) assess the association between outlier loci and climate. Using two primer sets, we genotyped 378 plants at 384 polymorphic AFLP loci across regional ecotypes from central and eastern Kansas and Illinois. Neighbour-joining tree and PCoA revealed strong genetic differentiation between Kansas and Illinois ecotypes, which was better explained by IBE than IBD. We found high genetic variability within prairies (80%) and even fragmented Illinois prairies, surprisingly, contained high within-prairie genetic diversity (92%). Using Bayenv2, 14 top-ranked outlier loci among ecotypes were associated with temperature and precipitation variables. Six of seven BayeScanFST outliers were in common with Bayenv2 outliers. High genetic diversity may enable big bluestem populations to better withstand changing climates; however, population divergence supports the use of local ecotypes in grassland restoration. Knowledge of genetic variation in this ecological dominant and other grassland species will be critical to understanding grassland response and restoration challenges in the face of a changing climate.

Cold adaptation shapes the robustness of metabolic networks in Drosophila melanogaster.

When ectotherms are exposed to low temperatures, they enter a cold-induced coma (chill coma) that prevents resource acquisition, mating, oviposition, and escape from predation. There is substantial variation in time taken to recover from chill coma both within and among species, and this variation is correlated with habitat temperatures such that insects from cold environments recover more quickly. This suggests an adaptive response, but the mechanisms underlying variation in recovery times are unknown, making it difficult to decisively test adaptive hypotheses. We use replicated lines of Drosophila melanogaster selected in the laboratory for fast (hardy) or slow (susceptible) chill-coma recovery times to investigate modifications to metabolic profiles associated with cold adaptation. We measured metabolite concentrations of flies before, during, and after cold exposure using nuclear magnetic resonance (NMR) spectroscopy to test the hypotheses that hardy flies maintain metabolic homeostasis better during cold exposure and recovery, and that their metabolic networks are more robust to cold-induced perturbations. The metabolites of cold-hardy flies were less cold responsive and their metabolic networks during cold exposure were more robust, supporting our hypotheses. Metabolites involved in membrane lipid synthesis, tryptophan metabolism, oxidative stress, energy balance, and proline metabolism were altered by selection on cold tolerance. We discuss the potential significance of these alterations.

Heritable variation in garter snake color patterns in postglacial populations.

Global climate change is expected to trigger northward shifts in the ranges of natural populations of plants and animals, with subsequent effects on intraspecific genetic diversity. Investigating how genetic diversity is patterned among populations that arose following the last Ice Age is a promising method for understanding the potential future effects of climate change. Theoretical and empirical work has suggested that overall genetic diversity can decrease in colonial populations following rapid expansion into postglacial landscapes, with potential negative effects on the ability of populations to adapt to new environmental regimes. The crucial measure of this genetic variation and a population's overall adaptability is the heritable variation in phenotypic traits, as it is this variation that mediates the rate and direction of a population's multigenerational response to selection. Using two large full-sib quantitative genetic studies (N(Manitoba) = 144; N(South Dakota) = 653) and a smaller phenotypic analysis from Kansas (N(Kansas) = 44), we compared mean levels of pigmentation, genetic variation and heritability in three pigmentation traits among populations of the common garter snake, Thamnophis sirtalis, along a north-south gradient, including a postglacial northern population and a putative southern refuge population. Counter to our expectations, we found that genetic variance and heritability for the three pigmentation traits were the same or higher in the postglacial population than in the southern population.

Natural selection drives clinal life history patterns in the perennial sunflower species, Helianthus maximiliani.

In plants, ecologically important life history traits often display clinal patterns of population divergence. Such patterns can provide strong evidence for spatially varying selection across environmental gradients but also may result from nonselective processes, such as genetic drift, population bottlenecks and spatially restricted gene flow. Comparison of population differentiation in quantitative traits (measured as Q(ST) ) with neutral molecular markers (measured as F(ST) ) provides a useful tool for understanding the relative importance of adaptive and nonadaptive processes in the formation and maintenance of clinal variation. Here, we demonstrate the existence of geographic variation in key life history traits in the diploid perennial sunflower species Helianthus maximiliani across a broad latitudinal transect in North America. Strong population differentiation was found for days to flowering, growth rate and multiple size-related traits. Differentiation in these traits greatly exceeds neutral predictions, as determined both by partial Mantel tests and by comparisons of global Q(ST) values with theoretical F(ST) distributions. These findings indicate that clinal variation in these life history traits likely results from local adaptation driven by spatially heterogeneous environments.